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We uncover a nontrivial signature of the hierarchical structure of quasi-degenerate random di- 
rected polymers (RDPs) at zero temperature in 1 + 1 dimensional lattices. Using a cylindrical 
geometry with circumference 8 < W < 512, we study the differences in configurations taken by 
RDPs forced to pass through points displaced successively by one unit lattice mesh. The transi- 
tion between two successive configurations (interpreted as an avalanche) defines an area 5*. The 
distribution of moderatly sized avalanches is found to be a power- law P(S) dS ~ dS. Using 

a hierarchical formulation based on the length scales Wz (transverse excursion) and the distance 
Wz" between quasi-degenerate ground states (with < a < 1), we determine fj, = |, in excellent 
agreement with numerical simulations by a transfer matrix method. This power-law is valid up to a 
maximum size 5*5 ~ . There is another population of avalanches which, for characteristic sizes 

3 

beyond Ss, obeys P(S)dS ~ exp( — (5/Ss ) 3 ) dS also confirmed numerically. The first population 
corresponds to almost degenerate ground states, providing a direct evidence of "weak replica sym- 
metry breaking", while the second population is associated with different optimal states separated 

2 

by the typical fluctuation W s of a, single RDP. 
PACS numbers: 02.50.Ey, 05.70.Ln, 64.60.Ht 



I. INTRODUCTION 

Self-organized criticality (SOC) jl| describes out-of- 
equilibrium extended systems driven infinitely slowly 
which respond intermittently with avalanches or bursts of 
sizes distributed according to power-law distributions. A 
close relationship between critical phase transitions and 
a class of SOC systems |||| has been pointed out. Mem- 
ber systems of this SOC class operate exactly at the crit- 
ical value of an underlying critical point. A necessary 
condition for this to occur is that the order parameter 
(often akin to a flux) of a dynamical critical transition 
be driven infinitely slowly, thus forcing the control pa- 
rameter to readjust itself dynamically around its critical 
value ||. 

Motivated by this correspondence, we introduce a new 
SOC model. It can be described as an equilibrium de- 
pinning problem wherein a certain type of avalanche sep- 
arates local equilibrium states. The succession of equi- 
librium state transitions found in our model resembles 
the behaviour of abelian sandpiles Qj. In the latter, 
each avalanche can be shown to connect two different 
microscopic metastable states. Furthermore, its critical 
state is then characterized by the complete set of these 
avalanche-connected metastable states. Whereas the set 
of coexisting metastable stables is created by the thresh- 
old rules of the sandpile automata, the many coexisting 
local equilibrium states appearing in our model emerge 
from an optimal (i.e. minimum energy) configuration 
in a quenched random landscape. This disorder induces 



the coexistence of an extremely large number of almost 
equivalent configurations. The resulting closeness in en- 
ergy space leads to a large spread in configuration space. 
This turns out to produce a power-law distribution for 
the interconnecting avalanches. 

Thus, a common property of SOC systems is that they 
are characterized by a large set of almost equivalent and 
degenerate states. This set can be generated by dynamic 
automata rules, disorder, frustration or other mecha- 
nisms. In addition to the introduction of a new class 
of SOC models, our results provide further evidence for 
the hierarchical structure of sets of random directed poly- 
mers (RDPs). 

Our results bear an apparent strong similarity to those 
previously obtained for pinned charged density waves || , 
driven interfaces in random media , and elastic mani- 
folds on disordered substrates 0. However, the connec- 
tion between the dynamic critical phenomena obtained 
from a constant driving force F at the depinning thresh- 
old F c and the nearly critical behavior obtained by a 
small constant velocity drive is based on an argument 
relating the critical behavior as F — > F+ and F — ► F~. 
This predicts a vanishing exponent for the avalanche 
distribution in our 1+1 dimensional case, which seem- 
ingly contradicts our result. The discrepancy stems from 
the fact that we do not describe the same regime; the van- 
ishing exponent refers to the existence of large avalanches 
of sizes controlled by the system size (or the correlation 
length when off-criticality applies). This corresponds to 
the second of two identified avalanche regimes of our 
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model. In contrast, the present work reveals the existence 
of a sub-dominant power-law distribution of avalanches 
stemming from the hierarchy of almost equivalent degen- 
erate states. These states do not, however, contribute to 
the large scale behavior and have thus been overlooked 
in previous work. 




The model is defined in the next section, while in 
Sec. Ill we derive our theoretical predictions for the dis- 
tribution of avalanche sizes. These are compared with 
extensive numerical simulations in Sec. Our conclu- 
sions are found in Sec. [v| 
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FIG. 1. A typical set of optimal configurations for a RDP of length W = 4096 and for < y < 1200. a) global system (grey 
framed boxes outlines regions of succeeding plots such that the horizontal and vertical extensions of these boxes follow Eqs. (H) 
and ([]) with a » 0.9); b) magnification of the largest box in a); c) magnification of the largest box in b); d) magnification of 
the box in c). 



II. DEFINITION OF THE MODEL 

Consider a RDP on a square lattice oriented at 45° 
with respect to the x axis and such that each bond carries 
a random number, interpreted as an energy. An arbitrary 
directed path (a condition of no backwards turn) along 
the x-direction and of length W (in this direction) corre- 
sponds to the configuration of a RDP of W bonds. In the 
zero temperature version we study here, the equilibrium 
polymer configuration is the particular directed path on 
this lattice which (in the presence of given boundary con- 
ditions) minimizes the sum of the W bond energies along 
it. This simple model, with its much varied behavior, has 



become a valuable tool in the study of self-similar surface 
growths H , interface fluctuations and depinning 0] , the 
random stirred Burgers equation in fluid dynamics |ic| l 
and the physics of spin glasses |ll]] . 

Let us apply a field h that exerts a force on one of the 
vertical endpoint position y(W) of the polymer. This 
field adds a term — hy(W) to the configurational energy 
of the polymer given by the sum of random bond energies 
along it. It is similar to a transverse electric field acting 
on the charged head of the polymer. If the other poly- 
mer extremity is free, the minimum energy is obtained 
by letting y(W) go to infinity as the external field term 
— hy(W) diverges to — oo. This energy always dominates 
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the configuration energy for any reasonable distribution 
of random bond energies. A depinning transition thus 
occurs for the value h — + of the control parameter h. 
Mezard [O] has shown that holding the other endpoint 
fixed results (in the small field limit) in extremely jerky 
displacement of the charged head as a function of the field 
strength h. The position of the charged head is stationary 
for large ranges of applied field values and then changes 
suddenly. At the field values where these transitions (or 
avalanches) occur, the susceptibility attains large values. 
These susceptibility bursts are reportedly distributed ac- 
cording to a power-law ]T^] . This avalanche response has 
been attributed to a "spin glass phase" , with several 
valleys of similar energy. It is important to realize that 
this avalanche behavior is not SOC as the driving is non- 
stationary; nothing occurs when the field stays constant 
and increasing the field will lead ultimately to the situa- 
tion where the RDP is blocked in a fully extended con- 
figuration along the first quadrant bisectrix. This regime 
is similar to a mode of operation with a slow sweeping of 
a control parameter 111}] . 



polymer as wrapping itself around a cylinder of circum- 
ference W. The polymer is allowed to equilibrate, i.e. 
take the spatial configuration of minimum total bond en- 
ergy. Let us now shift the vertical position of the fixed 
end points from y to y + 1 (where the lattice mesh is 
taken as unity). The polymer is again allowed to equi- 
librate to the spatial configuration of minimum energy. 
We continue in this fashion in an iterative process, which 
amounts to controlling the average vertical velocity of the 
polymer to a value so small that the time scale to move 
over a unit mesh is much larger than any relaxation times. 
This guarantees that the polymer always finds the spatial 
configuration of minimum bond energy. Note that the y 
position of the polymer end points therefore functions as 
a clock, since no other relevant time scales are present. 

Figure ^ shows a typical set of optimal configurations 
for a RDP of length W = 4096 and for < y < L = 1200. 
The magnifications illustrate the self-affine structure of 
the RDPs and the self-similar hierarchical pattern of the 
local branching structure. 




FIG. 2. Schematic representation of optimal RDPs fixed 
at their two endpoints. An avalanche is defined by the area 5* 
spanned by the transition from the optimal configuration at 
y to y + 1, i.e. S is the area interior to the perimeter formed 
by the union of the two optimal RDP configurations at y 
and y+1 and the two vertical segments [(0, y); (0, y + 1)] and 
[(W, y); (W,y + 1)]. The successive avalanches are represented 
in different grey scales. 

The correspondence between depinning transitions and 
SOC models || suggests, with the previous results |l2| , 
the following variant of the problem. Instead of applying 
a field (control parameter), we set the depinning velocity 
(order parameter) to an infinitesimal value Jl4| ]. This is 
accomplished by initially fixing the two ends of the poly- 
mer at (x\ = 0, j/i = y) and {x2 = W, y<i = y). Since the 
two ordinates y\ = yi ~ y are equal, we could consider 
the case where only one endpoint is fixed while keeping 
the other one free, therefore making this situation cor- 
respond to a polymer on average twice as long but with 
both endpoints fixed. Alternatively, we may consider the 



III. THEORETICAL PREDICTION OF THE 
AVALANCHE SIZE DISTRIBUTION 

In this model, an avalanche at y is simply the transi- 
tion from the optimal configuration of a RDP with end 
points fixed at yx = y<i = U to the optimal configuration 
where the end points are now at yx = y% = y + 1 (as 
shown in Fig. g). We define the size of an avalanche by 
the area S spanned by the transition from the optimal 
configuration at y to the one at y + 1, i.e. S is the area 
interior to the perimeter formed by the union of the two 
optimal RDP configurations at y and y+1 and the two 
vertical segments [(0,y); (0, y+1)] and [(VF,y); (IF,y+l)] 
(see Fig. ||) . What do we know about the distribution of 
these avalanche sizes? 

Clearly, the structure of the ensemble of the optimal 
RDPs (for all possible end point y locations) uniquely 
determines the avalanches. For a RDP of length W, it is 
known that the typical transverse excursion Y in 1 + 1 di- 
mensions scales as Y ~ W v with v — | (see jll| and ref- 
erences therein). We thus expect that there exists a class 
of RDP transitions with vertical lengths Y of at least the 
order of this typical transverse excursion. The area S 
spanned by such a transition is therefore proportional to 
S ~ WY ~ W 3 (i.e. the characteristic avalanche size). 
The distribution of Y is known to behave asymptotically 
as P(Y) ~ exp(-(^r) 3 ) 111. Subtituting for S ~ WY 

in P(Y) give us 



P(S) ~ exp(- 



S 



(1) 



for S at least of the order of W 3 . This constitutes our 
first prediction. Its validity will be tested numerically in 
SecllV. 
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We now derive the distribution of avalanches in the 
large W limit for 5* smaller than . First notice that 
the sequence of optimal paths with fan shaped families 
of end points strongly resembles the ranking of paths by 
Zhang |lq] . He found that the difference Y between the 
endpoints of these optimal paths scale with path length 
X as Y ~ I" s with v a = |. This property is important 
for the understanding of our results as it suggests a hier- 
archical structure. We thus briefly recall its derivation. 

The Bethe ansatz with the replica trick Jl?]] provides 
a solution of the RDP problem in 1 + 1 dimension. This 
shows that the RDP problem is equivalent to solving a 
problem of n bosons in one spatial dimensions interact- 
ing with an attractive delta function potential. In this 
framework, imposing conditions on the endpoints of the 
RDP implies that the Bethe ansatz wave function must 
incorporate the motion of the center of mass of the n 
bosons: 

1 " 

*~l/exp(^|a; Q -^| + — ■ (2) 

ct,/3 a=l 

The term x 2 a represents the kinetic and 

3 \ Xa ~ x 0\ * ne potential energy. Since x a ~ wi , 
the kinetic energy ~ W^. In the Bethe ansatz wave 
function, the potential energy must be comparable to 
the kinetic energy, thus \x a — xp\ ~ , confirming that 
v s = |- This scaling describes the distance between de- 
generate ground states with so-called "weak replica sym- 
metry breaking" |l7| . Technically there is a replica sym- 
metry breaking but the distance between the degenerate 
states becomes negligible compared to their intrinsic fluc- 
tuations in the thermodynamic limit. 

We generalize the above observations to infer two 

2 2 

transverse length scales W$ and Wz a , where < a < 1 
(the case a = 1 is addressed separetely below), to de- 
scribe the hierarchical structure of RDP configurations 
as exemplified by Fig. Our results turn out to be inde- 
pendent of the choice of a. Intuitively, a family of width 
Wz a consists of families of width Ws a each of which 
consists of families of smaller width and so on (down to 
the elemental scale of the mesh) . The width, number and 
other properties of these embedded sets of families can 
be obtained from the two length scales Wz and W^ a 
using only dimension conservation and self-similarity ar- 
guments. 

1. The highest order family, that we call of order 1, 
corresponds to all the locally optimal paths that are 
within a distance of order of a best path. The 
vertical width of this family of order 1 is wi oc \Vs . 
This family is composed of locally optimal paths 
that join after a distance l± oc W, obtained by 

2 2 

the condition that If oc W ^ (this condition will 
become nontrivial at lower levels of the hierar- 
chy). The generic area covered by this family is 
Si oc l\ wi ~ Wz. This is also the typical size 



of the largest possible avalanche as defined above 
and corresponds to a transition between members 
of this family of order 1 . 

2. Within this family of order 1, we define N 2 fami- 
lies of order 2, each of which have a characteristic 

2 Q 2 

width w 2 oc 11 ~ Wz a . It is at this point that we 
have used the second length scale introduced by the 
quasi-degenerate ground states. From the conser- 
vation of (vertical) width, we have by construction, 

N 2 w 2 = wi , (3) 

leading to N 2 oc wi^~ a \ A family of order 2 
is by itself composed of locally optimal paths that 
join after a distance l 2 oc W a , obtained by the self- 
consistent condition that 

2 

w 2 oc Z| . (4) 

As a consequence, the generic area, i.e. the largest 
possible avalanche, covered by this family of or- 
der 2 (intra member transitions) is S 2 oc l 2 w 2 ~ 
W a x Wi a = Wi a . 

n. We infer that the relevant quantities of the n-th or- 
der family only depends on the associated ones in 
the family of order n— 1. This leads us to a recursive 
scheme for the calculation of the above introduced 
entities. In what follows we will formally define the 
simplest version of the iterative system of equations 
and state its solutions. 

Within each of the families of order n, we define N n+ i 
families of order n + 1, each of which have a characteristic 
(vertical) width ui n +i. From the conservation of width, 
we have by construction, 

N n+ iw n+ x = w n . (5) 

The characteristic width w n+ i relates the generic dis- 
tance l n after which locally optimal paths (within a fam- 
ily of order n) typically join. It obeys 

2 

w n+1 = aln a ■ (6) 
But the self-consistency condition relates w n+ i to l n +\ 

wn+i = Blf l+1 . (7) 
We are thus led to the direct recursion 

ln+l = All ■ (8) 

The typical area covered by an avalanche among the fam- 
ilies of order n + is 

S n +i — Cl n +iw n +i ■ (9) 
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Here A, a, B, and C are (real valued) constants. Since 
li oc W, we also have an initial condition for the recur- 
sion. This is generalized as 



with an exponent 



h = f(w) . 



(10) 



Finally, we have for the total number of families A/"„+i 
up to and including order n + 1 , 



M n +1 = N n+1 N n 

M = i • 

We find from Eqs. (§) and @ 



.4 



and from Eq. (fjj) that 



for n > 1 



Together with Eq. (|5|) we then get 

N n+1 

and from Eq. (^| 

S n = CB- 



rf(w)y {1 - a)0!T ' 



The latter expression is conveniently turned into 



5 



S n 



fJWY ) 



(11) 
(12) 



(13) 



(14) 



(15) 



(16) 



(17) 



For the cumulative number of families to order n, i.e. 
Eqs. Q and @, we get with Eq. @ 



(18) 



which with Eq. (17) results in a direct 5„ and VF depen- 
dence, 



■A/"„ = (CB)^/(W)3/^ 



(19) 



This reasoning, based on the hierarchical model, gives 
us the number of avalanches of specific sizes. To get the 
probability density distribution, we have to divide this 
number by the interval width from S n to S n +i which 
is simply proportional to S n up to a correction of order 
as seen from Eq. (|l6|). Gathering all the pieces 
and assuming that f(W) oc W leads us to the following 
prediction for the distribution of avalanche sizes 



P{S) dS oc dS , 



(20) 



M= 5 



(21) 



The power law in Eq. (|20|) describes the distribution of 

avalanche sizes 1 < 5 < Ss (we define Ss oc W^) in 

3 3 
the limit W — * oo. This upper scale Ss corresponds to 

the maximum typical sizes of the avalanches of order 1 

in the hierarchy. Notice that the prediction of Eq. ( |2l| ) 

is independent of the value < a < 1 and is thus robust 

with respect to the detailed structure of the hierarchy. 

A similar hierarchical structure can also be constructed 
for a = 1. In this case, it is postulated that u>„+i = 
w n /\, where A > 1 is the constant reduction factor from 
one level of the hierarchy to the next. While keep- 
ing Eq. (||), this leads to l n+ \ — l n /X^ and thus to 
Sn+i = Sn/X 1 using Eq. ([)]). The total number of fami- 
lies of order n is now simply proportional to A™. Solving 
as a function of S n , we retrieve exactly expression ( |20| ) 
for the distribution of avalanche sizes. 

This derivation is simpler because the hierarchical 
structure is exactly self-similar, with the same scaling ra- 
tio A throughout. This is in contrast to the case a < 1 for 
which Aocw 3 " ( 1-Q ) decreases with increasing family 
order. This derivation for a = 1 and A > 1 clarifies 
the origin of the exponent fi = 2/5 stemming simply 
from l//i = 1 + 3/2, i.e. from the fundamental self-aflme 
structure of the RDP with transverse excursion exponent 
2/3. 

We thus stress that the prediction of Eqs. (^0|) and 
( pl| ) is very general and independent of the specific hier- 
archical structure of the sub-dominant quasi-degenerate 
ground states. Note that the power-law distribution 
given by Eq. ( ^p| ) for the spanned surfaces is associated 
with two other power-laws, namely that for the distri- 
bution of typical transverse deviations w and that for 
the distribution of typical longitudinal deviations I. This 
stems from S ~ wl and w ~ Z 2 / 3 , leading to 



P(w) dw 



dw 



w 



l+i 



and P{1) dl 



dl 



11+2/3 



(22) 



As mentioned in Sec. Q, our finding seem to be in dis- 
agreement with the predicted value (equal to zero) for 
the avalanche size distribution in 1 + 1 dimensions . 
However, the avalanche regime of interest to us is differ- 
ent from that previously investigated. Our regim consists 
of small and intermediate avalanches of sizes up to Ss , 
whereas the regime which contains avalanches larger than 
Ss yields a size distribution with a vanishing exponent. 
The present work proposes a sub-dominant power-law 
distribution of avalanches which originates in a hierarchi- 
cal ordering of the almost equivalent degenerate states. 
Previous work has addressed the tail end of the avalanche 
distribution and has therfore not been attentive to the 
presence of these states. 
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IV. NUMERICAL TESTS 



The distribution of avalanche sizes S has been deter- 
mined numerically by using a now standard transfer ma- 
trix method Jlq ] relying on the chain property applying to 
the energy e(aq, yi\ £2,2/2) of a RDP going from (xi,yi) 
to (2^2,2/2): 

e(xi,?/i;a:2,?/2) = mm[e(x 1 ,y 1 ;x' ,y') + e(x' , y'; x 2> 2/2)] ■ 



(23) 



Figure || shows the distribution of avalanche sizes ob- 
tained numerically for system widths from W — 8 to 
512. For each width, we have calculated the RDP config- 
urations and the corresponding avalanche areas for sys- 
tem lengths 3 x 10 6 < L < 2 x 10 8 . These very long 
systems provide reliable statistical estimates. In Fig. |3j, 
the existence of a power-law for the distribution P(S) is 
quite apparent. The size interval over which the power- 
law holds increases as Ss ~ W 5 ^ 3 . Another feature of 
Fig. [| is the clear evidence of a characteristic avalanche 
size, corresponding to the bump of the distribution in 
the region of large avalanche sizes. The location of these 
bumps also scales as 5 s ~ W 5 ^ 3 . 
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FIG. 4. P(S) as a function of the rescaled variable S/W 5/3 
for W — 8 to 512 in a log-log plot. 
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tained numerically for system widths from W — 8 to 512 in 
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Finite size effects turn out to be very important in 
this problem and a careful finite size scaling analysis is 
appropriate. We approach this as follows. 
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For each system size, we determine the exponent fJ.(W) 
which best fits the numerical distribution. To demon- 



G 



strate the quality of the fit, we replot Fig. || by showing 
in Fig. H the function P(S) as a function of the rescaled 
variable S/Wi. For each size W, a different exponent 
n(W) is found. The dependence of (J,(W) as a function 
of W~ 3 is shown in Fig. ||. We find a very good fit ("least 
squares") with the finite size equation 



fj,(W) = Hoo - — - , 



(24) 



where c = 2.90 is a constant and /loc = 0.40. This is in 
excellent agreement with the prediction 2/5. An appar- 
ent power-law dependence of the exponent fi(W) on W 
as in Eq. (24) could result from fluctuations in the value 
of a within each level and accross the different levels of 
the hierarchy. 
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FIG. 6. P{s) for the different system sizes as a function 
of the rescaled variable S 3 /W 5 in a semilog plot. 

Finally, in Fig. || we represent P{s) for the different sys- 
tem sizes as a function of the rescaled variable S 3 /W 5 . 
This choice of variables is intended to test the prediction 
of Eq. |l]) . We observe a rather convincing tendency for 
the plot to converge to a straight line for large system 
sizes. 



V. CONCLUSION 

We have proposed a novel quasi-statically driven model 
that exhibits responses similar to those of SOC models. 
This model of a succession of optimal RDP configurations 
exhibits a power-law distribution of the area swept by a 
polymer between two successive optimal configurations 
(defined as an avalanche). 

Based on the existence of two fundamental scales W§ 
and W^ a (0 < a < 1) for the transverse fluctuations 
of a RDP of length W, we have constructed a hierarchi- 
cal representation of the set of quasi-degenerate optimal 
configurations. This hierarchy allow us to calculate ex- 
plicitely the exponent of the avalanche distribution. 



Our numerical analysis confirm the existence of two 
distinct populations of avalanches. One of these popula- 
tions consists of "small" avalanches that are distributed 
according to a power law with an upper cutoff controlled 
by the typical transverse length scale W$. The other 
population comprises the "large" avalanches beyond this 
typical transverse excursion Wz . 
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